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I. INTRODUCTION 


Lately, the properties of strongly interacting matter 
in the presence of external magnetic fields has become 
a very active research field. One of the driving moti¬ 
vations behind this interest is the possibility to study 
experimentally such properties in peripheral collisions of 
heavy nuclei at high energy. In addition, recent lattice 
QCD (LQCD) results show that the critical tempera¬ 
ture for deconfinement/chiral symmetry restoration de¬ 
creases with increasing field strength [l|. This behavior 
is dubbed inverse magnetic catalysis, and it reveals an 
unexpected, non-trivial phenomenon: in a thermal envi¬ 


ronment, near the transition temperature, the presence 
of a magnetic field acting on strongly interacting matter 
hinders the formation of a quark-anti-quark condensate. 
LQCD calculations @ show that the quark condensate 
does increase with increasing magnetic field at low tem¬ 
peratures. This behavior corresponds to magnetic catal¬ 
ysis. However, as the temperature increases approaching 
the crossover region T ~ 150 MeV, the quark conden¬ 
sate reaches a maximum value smaller than at T = 0 
(for the same field strength). Subsequently, the conden¬ 
sate decreases as a function of the field strength. Fi¬ 
nally, for temperatures above the cross-over values the 
condensate decreases monotonically as a function of the 
magnetic field. Some of the possible scenarios aim- 
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ing to understand this behavior include (i) invoking a 
fermion paramagnetic contribution to the pressure with 
a sufficiently large magnetization [3], (ii) the competi¬ 
tion between the valence and sea contributions at the 
phase transition Q produced by a back reaction of the 
Polyakov loop, which depends on the magnetic field 0, 
(iii) magnetic inhibition due to neutral meson fluctua¬ 
tions in a strong magnetic field Q, (iv) accounting for 
non-perturbative effects by means of Schwinger-Dyson 
Equations and renormalization group analyses 0-1, (v) 
a decreasing magnetic field and temperature dependent 
coupling with llOlIl2| and without [lj,[lj| plasma screen¬ 
ing effects, (vi) the proper account of the gluon contri¬ 
bution during the phase transition 0, and (vii) quark 
anti-screening due to the anomalous magnetic moment 
of quarks in strong fltSj and weak 0 fields. On the 
other hand, this behavior is not obtained in mean field 
approaches describing the thermal environment [HH23j], 
nor when calculations beyond mean field do not include 
magnetic effects on the coupling constants [24|. For re¬ 
cent reviews see Refs. [251 1261]. 

Given the dual nature of the QCD phase transition, 
a pertinent question is to what extent inverse magnetic 
catalysis is due to the mechanisms of either chiral sym¬ 
metry restoration and/or of deconfinement. One way to 
address this question is to find a relation between decon¬ 
finement and chiral symmetry restoration parameters as 
a function of the magnetic field. Since the transition hap¬ 
pens for temperatures in the realm of non-perturbative 
phenomena, the relation searched for needs to carry non- 
perturbative information. An extensively used tool in 
the context of effective models at finite temperature and 
zero (HI and finite (HI magnetic field is the Polyakov 
loop [ 23 . When coupled to quark degrees of freedom 
this loop sheds light on how chiral symmetry and decon¬ 
finement behave during the QCD transition as a function 
of the field intensity. Another non-perturbative tool that 
does not rely on effective models is that of QCD Finite 
Energy Sum Rules (FESR). This approach has been suc¬ 
cessfully applied both at zero (29j and at finite temper¬ 
ature [30( to understand hadronic properties. Of partic¬ 
ular mention are (i) the prediction, in this framework, 
of the survival of charmonium and bottonium above the 
critical temperature [H| . (ii) the temperature behavior of 
the hadronic width of charmonium and bottonium (31j 
in qualitative agreement with bottonium results from 
LQCD 0, and (iii) the description of the di-muon spec¬ 
trum in heavy-ion collisions [331 ]. in excellent agreement 
with data |3J] in the region of the rho-meson peak. A 
key parameter that emerges from this analysis signalling 
quark-gluon deconfinement is the squared energy thresh¬ 
old, soj above which the hadronic spectral function is 
well approximated by perturbative QCD (pQCD). An in¬ 
teresting relation between so and the quark condensate 
(qq ), whereby the former is proportional to the latter has 
been found, in the absence of a magnetic field and at fi¬ 
nite temperature, in [ 35 J, and at finite temperature and 



FIG. 1: Cauchy integration contour in the complex squared 
energy s-plane used to obtain QCD FESR. The radius of the 
circle, so, is the threshold for pQCD. 


density in f3Gl] . 

In this paper we use FESR in the axial-vector chan¬ 
nel, and in the presence of an external magnetic field, 
to explore the relation between (i) the deconfinement 
and chiral symmetry restoration parameters, So and ( qq ), 
and (ii) obtain the behavior of the gluon condensate as 
a function of the magnetic field intensity at zero temper¬ 
ature. Other formulations of the QCD sum rules, also 
at T = 0, have been used to explore the dependence of 
heavy-quark meson-masses and mixing as a function of 
the magnetic field strength [37J- The paper is organized 
as follows: In Sec. im we set up the FESR. We show that 
since the pQCD contribution to the current correlator re¬ 
ceives corrections that can be expressed as powers of the 
magnetic field strength divided by powers of the squared 
energy s, there appear additional terms that contribute 
to higher order FESR. In Sec. IIIII we explicitly compute 
the pQCD corrections to the current correlator and solve 
the first two FESR to find the dependence of So and the 
gluon condensate on the magnetic field strength. We 
show that for magnetic field strengths eB smaller than 
so, the former follows the magnetic field dependence of 
the quark condensate, which we parametrize from LQCD 
results 0 ]. We also show that the magnetic field depen¬ 
dence of the gluon condensate receives non-trivial correc¬ 
tions from the pQCD sector and that overall it is a mono- 
tonically increasing function of the field strength. We 
summarize and conclude in Sec. m leaving for the ap¬ 
pendices the explicit computation of the imaginary part 
of the hadronic contribution and the vanishing of the first 
order correction in the field strength of the pQCD con¬ 
tribution to the axial-vector current correlator. 
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II. FINITE ENERGY QCD SUM RULES IN THE 
PRESENCE OF A MAGNETIC FIELD 


The charged axial-vector current correlator in the ab¬ 
sence of a magnetic field and at T = 0 can be written 
as 


n^(g 2 ) = i J d A x e iqx <0\T(A^x), 4(0))|0 > 

= (~g^q 2 + QnQu) n A (q 2 ) + q^q v n 0 ((? 2 ), 

(i) 

where A^{x) =: d(x) r y IJj ^u{x) : is the (charged) axial- 
vector current, q^ is the four-momentum carried by the 
current, with s = q 2 > 0 the squared energy. The func¬ 
tions II J 4 i o(^ f2 ) are free of kinematical singularities, an 
important property needed in writing dispersion relations 
and sum rules. Concentrating on e.g. IIo(q 2 ) and invok¬ 
ing the Operator Product Expansion (OPE) of current 
correlators at short distances beyond perturbation the¬ 
ory, one of the two pillars of the QCD sum rule method, 
one has 

n o(Q 2 )lQCD = c 0 i + ^ (Q 2 Af(^ 2 )) , (2) 

N=1 

where Q 2 = -q 2 , ( 6 2 n{ g 2 )) = (0|O 2 iv(tt 2 )|0), and y 2 
is a renormalization scale. The Wilson coefficients C/v 
depend on the Lorentz indexes and quantum numbers of 
the currents, and on the local gauge invariant operators 
O 2 N built from the quark and gluon fields in the QCD 
Lagrangian. These operators are ordered by increasing 
dimensionality and the Wilson coefficients are calculable 
in pQCD. The unit operator above has dimension d = 0 
and CqI stands for the purely perturbative contribution 
normalized according to 

C 0 i=±\n(^j[l + O(a s (s))} , (3) 

Since there are no dimension d = 2 operators built from 
the QCD fields, it is generally assumed that the OPE 
starts at dimension d = 4. This is fully confirmed by de¬ 
terminations of condensates from experimental data [H[ ■ 
The dimension d = 4 in the chiral limit is proportional 
to the renormalization group invariant gluon condensate 

C 4 (6 4 ) = ^(a s G 2 ) (4) 

The second pillar of the QCD sum rule method is to 
consider an integration contour in the complex square 
energy plane, as in Fig. 1, and invoke Cauchy’s theorem 
assuming that QCD can be used on the circle of radius 
| so I, provided |so| is large enough (quark-hadron duality). 
On the real axis there is a discontinuity associated with 


the hadronic states entering the spectral function. Since 
there are no further singularities this leads to the FESR 

— -—: (f ds s Ar_1 IIg CD (s) = — f ds s Ar_1 ImIlQ AD (s),(5) 
27r * J C( |so |) ^ Jo 

with N > 1, and IIg CD (s) given by the OPE, Eq. ©. It 
will be shown later that in the presence of a magnetic 
field, and in the weak field limit eB < so, the Wilson 
coefficients acquire themselves a B-field dependence. In 
this work we shall compute the corrections to the FESR 
due to a weak magnetic field, which can be expressed as 
a series in powers of eB. Since the magnetic field carries 
dimension of energy squared, on dimensional grounds one 
finds the replacements 



-A 



£e<" 

n —1 


) (eBY 


C'2N 


sr r <">> {eB)m 

/ °2JV m 
m —0 


( 6 ) 


where C are dimensionless quantities that can be com¬ 
puted in pQCD at a given order in eB. Substituting 
Eqs. © and © into Eq. ©, one obtains 


JV-l 

-£(-l ) N - m C^N- m ^Br(0 2{N - m) ) = 

m —0 

- [ S dss N -HraIl^(s)-^ + d 0 N \eB) N . ( 7 ) 
n Jo -'V 

Note that in general the presence of the magnetic field 
mixes operators of different dimension in the FESR. For 
instance, the first two sum rules (N =1,2) become 

1 r s ° 

0 = ds Imn“ AD (s) - C 0 s 0 + C ( u\eB) , (8) 

n Jo 


I c 

-Cf ) ( 0 4 )+C«(eS)( 0 2 ) = - / 

II Jo 


ds slm n“ AD (s) 


- ^s 2 + C^(eB) 2 .( 9) 


In order to set up explicitly the relevant FESR we 
start by computing the hadronic contribution. The axial- 
vector current in the presence of a magnetic field can be 
interpolated by the charged pion current 

4 = “A- D M 7r+ = “ ieAp) 7T + , (10) 

where = 130.28(14) MeV [39j] is the pion decay con¬ 
stant, 7 r + the pion field, and A M = (B/ 2) (0, —y,x, 0) the 
vector potential in the symmetric gauge, which gives rise 
to a constant magnetic field along the z direction. There¬ 
fore, the axial-vector correlator in the hadronic sector can 
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be written as 

n ™{x,y) = < 0\T(A„(x) , 4(y))|0 > 

= if%(0\T [D^ir + (x)D*ir~(y)] |0> 

= if% D »(x) D t(y)Gn(x,y) 

= ie^'V) fl (d/dx?) 0 d/dy ") - y) 

= e«< s ^n“ D (x-!/), (11) 


where we have used the fact that the charged pion prop¬ 
agator Gtt(x, y ) in the presence of a magnetic field can be 
written as a product of a transnationally invariant piece 
G n (x — y) and a phase factor e* e$ ( x,!/ ). The phase factor 
does not depend on the integration path so that choosing 
a straight line trajectory it can be written as 

®(x,y)= [ AOdt ( 12 ) 

J X 


where £ = yt + x(l —t),t£ [0,1]. It is easy to show that 
the above phase factor can be gauged away by a suitable 
gauge transformation of the vector potential. Hence, we 
keep only the translational invariant part of the hadronic 
correlator whose Fourier transform is 


n rV) = i/X(9 2 ), (13) 

where G^(q 2 ) stands for the Fourier transform of Gn(x — 
y). Using Schwinger’s proper time method this quantity 
can be written as 

G n (q 2 ) = / - 7 - ir[ q ^-qlt a n(eBr)/eBr+ie] ( 14 ) 

J o cos (eBr) 

where m n = 0 as we consider the chiral limit. Hereafter 
we shall use the notation 


— ffjL 9fj.v 
sJL = diag(l, 0 , 0 , — 1 ) 

9»„ = diag(0,1,1,0). (15) 

Consequently 



a ■ b 

<3 

1 

Co' 

II 

(a 


= glaW 



= a 0 b 0 - a 3 b 3 

(a- 

b)x 

= 9^b v 



= ai&i + 0262 

J 

y in 

,g^ 

= 2 

9^ 

,9^ 

= - 2 . 


(16) 


Extending the integration in Eq. ED to the lower com¬ 
plex plane, it can be shown [401 ] that the charged pion 
propagator is expressed in terms of a sum over Landau 
levels 


G„(q 2 ) = 2lJ2 
1=0 


(-\) l U(2q\/eB)e~^l eB 
< 7 jj — (2 1 + l)eB 


( 17 ) 


As it is customary in problems where Lorentz invari¬ 
ance is lost and the correlators depend separately on the 
time and space parts of the four-momentum, to simplify 
the analysis hereafter we study the behavior of the cor¬ 
relator in the static limit where the space part of the 
four momentum vanishes. However, note that since the 
magnetic field separates space-time into longitudinal and 
transverse directions with respect to the direction of the 
field, q 0 comes together with q 3 to form gjj = q 3 — q 2 . 
It is therefore equally simple to set q\ = 0, perform the 
analysis in terms of qjj and then at the end, if needed, 
take <73 —>- 0. In this limit, Eq. ED becomes 


G^{q ) — 2i ^ 

i=o q \\ 


(- 1 Y 

(21 + l)eB ’ 


(18) 


Therefore the hadronic contribution to the correlator be¬ 
comes explicitly 

°° f-lV 

n 5“(4 = ,) = - 2 /ig J -l n L 5 . (is) 

As we show in Appendix A, the imaginary part of 
Eq. (fl9l) in the weak field limit, eB < sq, is given by 

Imn“ AD (s) = f 2 7rS(s - eB). (20) 

It is important to mention that Eq. (EUll holds provided 
so satisfies also the condition sq < 3 eB, as explicitly 
shown in Appendix A. This means that the weak field 
condition (eB < sq) has to be supplemented with a fur¬ 
ther restriction for Eq. (l20l) to remain valid. Recall that 
so represents the onset for the pQCD description for the 
axial-vector spectral density and that this quantity is a 
decreasing function of temperature [33j. For cold nuclear 
matter, as in the case of a neutron star, the condition 
eB < so < 3 eB may be difficult to meet, specially for a 
weak field where not only the situation eB < So, but even 
neB < so, (n > 1) can happen. However, for a heavy- 
ion collision, around the deconfinement/chiral symmetry 
restoration transition, when so has dropped off to small 
values, the weak field condition can also be made com¬ 
patible with eB < so < 3 eB. Hereafter we keep in mind 
this last observation as the working scenario, aiming to 
eventually incorporate thermal effects to describe the be¬ 
havior of a magnetized medium near the phase transition. 


Using Eq. (l20l) . the hadronic line integral in the FESR 
is given explicitly by 

- [ ds s JV_ 1 ImnQ AD (s) = fl (eB) N ~ x . (21) 

7r Jo 

Substituting Eq. (l2ll) into the QCD sum rules Eqs. 4HD- 
® gives 

0 = fl-Coso + C < d ) (eB) 

-C 4 (0 4 ) = fl(eB) - + Cf (eB) 2 . (22) 

In order to solve these equations, we now proceed to com¬ 
pute explicitly the coefficients and C^. 
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k-q 



FIG. 2: pQCD contribution to the axial-vector current corre¬ 
lator in the presence of a magnetic field. The thick internal 
lines represent the full quark propagators in the magnetic field 
background. 


III. SOLUTIONS FOR s 0 AND C 4 (0 4 ) 


To perform the perturbative calculation of the coeffi¬ 
cients Cq 1 ' 1 and C'g 2 ' 1 in the weak field limit we make use 
of the weak field expansion of the quark propagator in 
the presence of a constant magnetic field [41| , and in the 
chiral limit, up to order 0(B 2 ) 

iS B {k) = 


where e q is the absolute value of the quark’s charge. 
We emphasize that the weak field expansion of the the 
fermion propagator in the chiral limit is a well defined 
object. This happens because the weak field limit of the 
propagator can be thought of as a series representation 
in powers of eB of the full propagator independent of any 
relation between the strength of the magnetic field and 
the fermion mass. One can then use this series repre¬ 
sentation and take the chiral limit as much as one can 
do it for the usual fermion propagator in the absence 
of the field; depending on the process considered, this 
limit could lead to infrared divergences but it does not 
preclude the validity of taking m —> 0 inside the propa¬ 
gator. In the present context, the field can be considered 
as weak when compared to so, which is the only other 
energy (squared) present in the problem, in such a way 
that the expansion in Eq. © converges. 

The pQCD contribution to the axial-vector current 
correlator in the presence of a magnetic field is depicted 
in Fig. [2j where we also define the kinematics. The thick 
internal lines represent the full quark propagators in the 
magnetic field background. To first order in e q B only 
one of the two quark propagators carries the magnetic 
effects. This is depicted in Fig. [3] where the wavy line 
starting from a cross represents the external magnetic 


. # , m 7172(7 •*) || 


'k 2 
2 i(e q B) 2 


fc 4 


*i(7-fc)|| 'k) 


k 8 


(23) 



FIG. 3: pQCD contribution to the axial-vector current corre¬ 
lator in the presence of a magnetic held to first order in e q B. 
The thick wavy lines ending in a cross represent the external 
magnetic held. 


field. The two diagrams in Fig. [3] that determine the 
coefficient C^\ vanish identically when contracted with 
the momenta carried by the axial-vector currents. This 
is due to a straightforward application of Furry’s the¬ 
orem, and to the fact that the vector and axial-vector 
correlators are chiral symmetric. However, we show ex¬ 
plicitly in Appendix B that each of the contributions in 
the diagrams, Fig. [3] vanishes. 

The first non-trivial magnetic contribution to the 
pQCD axial-vector current correlator is of order (e q B) 2 . 
The relevant diagrams are shown in Fig. [I] First, we 
compute the diagram where one magnetic field line is at¬ 
tached to each one of the two quark propagators. For 
these, we use Eq. C3l) to first order in e q B. We call this 
contribution n\}„\q 2 ), and its explicit expression is 

I^V) = - iN c {q u q d B 2 ) 

f d 4 k Tr [ 7^7172 [7 • (fc ~ g)||] 7*/7i72(7 • fc)||] 

J (27t) 4 (fc — g) 4 fc 4 

(24) 


Since according to Eq. ©, we are interested in the 
magnetic corrections to the coefficient of the longitudi¬ 
nal structure, no(< 7 2 ) we project H^\q 2 ) with q^q" and 
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FIG. 4: pQCD contribution to the axial-vector current cor¬ 
relator in the presence of a magnetic field to second order in 
e q B. The thick wavy lines ending in a cross represent the 
external magnetic field. 


define 

n^V) = ^n^fo 2 ). (25) 

Using 

7i72 [7 ' (k - g)||] = [7 ' (k - q)\\] 7i72 (26) 

together with 

7172 #7172 = 2(7 • q )_l - 4, (27) 

gives 

d 4 k 


I#' = iN c {q u q d B 2 ) 


TV 


, (2tt) 4 

(k~q)] || (—2(7 • q)± + 4){l ' ^)|| 


(k — q) 4 k 4 

The evaluation of the trace yields 
TV 


(28) 


• (k - g)]|| (-2(7 • q)± + 4 )(7 • fc)n 
= -4(5jj + «l)[(fc-5)-fc] | | 

+ 8(g- fc)|| [(fc — q) - 9] m ■ 


We now use the Feynman parametrization 


1 


= 3! / dx- 


x{x — 1 ) 


r> (30) 


C k-q) 4 k 4 ' J 0 [(fc - xq) 2 - x(x - l)q 2 Y 

and the change of variable 

k —>• l = k — xq , (31) 

to obtain 

n i 0 11) (q 2 ) = 4iN c {q u q d B 2 )3\ x(x - 1) J 1 4 

9 2 gjja;(a;-l)-(g| + gi)^ + 2(g-/)| , (32) 

where we have discarded terms with odd powers of l and 
defined A = x{x—l)q 2 . The integrals over l are computed 
by means of 


d d l 


1 


(2n) d [P - A] r 
d d l PP 
(2n) d [P - A] r 


= 1 


(-l) n T(n — d/2) / 1 \ n ~ d/2 

(47r) d / 2 r(n) \A/ 

.(-1 ) n ~ 4 T{n-d/ 2-1) 


X '1 


(47r)d/2 2 

n-d/ 2-1 


F(n) 


(33) 


with n = 4 and d = 4. Using Eq. (1531) in Eq. (1351) . and 
after integrating over x, we find 


1(11). N r- 


n r ; =-^(^ 2)L 


d\+q\ 


In the limit q\ —> 0 Eq. (1341) becomes 


n, 


( 11 ) ^0 n c 


-^(quqdS 2 )- 


(34) 


(35) 


In a similar fashion we compute the diagrams in Fig. [5] 
to second order in eB in the u-quark and in the d-quark 
propagator. Calling the longitudinal projections IIq 20 * 
and n(° 2 \ respectively, the result is 


nfV) = - 


n o 02) (d 2 ) = 


N c 
24t r 2 

gi->0 N c 
247r 2 

N r . 


247r 2 

gi->o N c 
> _ 24tt 2 


{quB) 2 
( q u B ) 2 
(q d B) 2 
(■ q d B ) 2 . 


(ql + qj) gjgj 

Q I | 


(d|| +9J.) j n qi 

o 1-^1 


(36) 


Adding all three contributions, and using the absolute 
values q u = 2/3 e, q d = 1/3 e, and N c = 3, we obtain 
the coefficient of the longitudinal structure of the axial- 
vector current correlator to second order in the magnetic 
field 


n o =-It= 


(eB) 2 
47r 2 


(29) 


17 

18 


(37) 
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FIG. 5: (Color on-line) Light-quark condensate normalized to 
its vacuum value as a function of the magnetic field strength 
eB in units of GeV 2 . The data points are from Ref. ;2j and the 
dotted line corresponds to the fit (qq) ( eB ) / (qq) = 1 + a(eB) + 
b(eB) 2 , with a = 0.85 GeV" 2 , b = 0.34 GeV" 4 . 


Using this result together with the first equation in 
Eq. ([ 6 ]), we obtain the Wilson coefficient of the pQCD 
contribution to second order in the magnetic field 

c » 2) =- GO b < 38 > 

The last ingredient needed to find sq and 64 ( 04 ) is 
the magnetic field dependence of This may be ob¬ 
tained from the Gell-Mann-Oakes-Renner (GMOR) rela¬ 
tion, which relates with the light-quark condensate, 
and the pion mass with the light-quark masses. The 
GMOR relation is a low energy theorem involving the 
pseudoscalar current correlator 

Ml 2 )=i [ d 4 xe^ <0\T(d^A^x)d’'Al(0))\0>, 

(39) 

where the divergence of the axial-vector current, 
d^A^(x), is proportional to the sum of the up- and down- 
quark masses. At zero momentum the GMOR relation 



0.4 0.5 0.6 0.7 0.8 0.9 1.0 

eB/so(0) 


FIG. 6: (Color on-line) Squared energy pQCD threshold, so, 
normalized to its eB = 0 value, as a function of the field 
strength eB. Note that the plotted range is consistent with 
the requirement eB < so < 3 eB. 


reads 


ml fl = ~ 2 (m u + m d ) (qq) , (40) 

In the standard chiral interpretation of this relation 
is related to the light quark condensate 

fl = ~2B(qq), (41) 

and the pion mass is related to the light-quark masses 

ml = ^(m u + m d ), (42) 

where B is a CHPT constant. In order to find the con¬ 
stant B one needs to use the physical pion and light-quark 
masses. In the hadronic sector the pion pole contribution 
to the pseudoscalar current correlator, ipbq 2 ) is given by 

V> 5 HAD (9 2 ) = */>"G.(g 2 ), (43) 

where G^(q 2 ) is the pion propagator. In the QCD sector 
the light-quark condensate contribution to ^ 5 ( 0 ) is given 
by 


V> 5 QCD ( 0 ) = -2 (m u + m d )(qq). (44) 

According to Eq. (1T81) [see also the discussion leading to 
Eq. (1561) for the normalization] the pion propagator for 
q 2 > 0 in the presence of a weak magnetic field can be 
approximated as 


G„(q 2 ) - 


— m 2 — eB 


(45) 


so that the GMOR relation becomes 

f 2 m\ 


eB 


— -2 (m u + m d )(qq). 


For weak fields, we use the expansion 


1 


eB 


m ^ m ^ 


mt 


(46) 


(47) 


Inserting Eq. (1471) into Eq. (l46l) . and neglecting a mag¬ 
netic field dependence of the pion and quark masses, the 
term introduces a correction to the GMOR relation of 
higher order than the one being considered in this work. 
Therefore, to a good approximation for small fields the 
magnetic field dependence of f 2 is determined by that of 
< qq > and the constant B from its vacuum value. We 
point out that the validity of the GMOR relation in the 
presence of a magnetic field has first been shown in the 
first of Refs. (Til. 


The light-quark condensate in the presence of the mag¬ 
netic field has been computed in Ref. jS}. We make use 
of this result, and parametrize the magnetic field depen¬ 
dence of the light-quark condensate with a quadratic fit 

(m)(eB)/(qq) = 1 + a ( eB ) + fo (e£) 2 , 


(48) 









0.4 0.5 0.6 0.7 0.8 0.9 1.0 

eB/s 0 (0) 


FIG. 7: (Color on-line) The gluon condensate, < 74 ( 04 ), nor¬ 
malized to its eB — 0 value, as a function of eB. Note 
that the plotted range is consistent with the requirement 
eB < so < 3 eB. 


where a = 0.85 GeV 2 , b = 0.34 GeV 4 and (eB) is given 
in GeV 2 . The data from Ref. 0 together with the fit are 
shown in Fig. [5] Using this information we finally write 
the explicit solutions for Sq and < 74 ( 04 ) from Eq. (1221) 


so 

o 4 (o 4 ) 


-8 tt B(qq)( eB ) 


-2(eB)B(qq) (eB) + 8tt (B (qq) {eB) ) 2 
^ 18 y 4 tt 2 ' 


(49) 


as a function of the magnetic field allows then to extract 
the magnetic field dependence of /„.. We should point 
out that since lattice QCD results are valid for physi¬ 
cal pion masses, the calculation in this work cannot be 
directly compared to the results in Refs. [HI which are 
performed using chiral perturbation theory in the limit 
where qB m 2 . The magnetic field dependence of So 
is proportional to the magnetic field dependence of the 
absolute value of the light-quark condensate. This means 
that the behavior of these two parameters as a function 
of the magnetic field is similar. Therefore the magnetic 
field both helps the formation of the condensate and acts 
against deconfinement. Next, we computed explicitly the 
magnetic corrections to the pQCD contribution which af¬ 
fect the behavior of the gluon condensate as a function 
of field strength. In the weak field limit the first correc¬ 
tion is quadratic in the field. The gluon condensate also 
grows as a function of the field strength which goes hand 
in hand with the behavior of the magnetic field, both as 
a catalyst of chiral symmetry breaking and confinement. 

The results obtained here should serve as a basis for 
studies at finite temperature in an external magnetic 
field, thus providing clues on the relation between chiral 
symmetry restoration and deconfinement at the transi¬ 
tion temperature. This research is in progress and will 
be reported elsewhere. 


The solutions for s 0 and for < 74 ( 04 ) as functions of eB 
are plotted in Figs. [S]and[7J respectively. Note that s 0 
is proportional to the absolute value of the light-quark 

condensate, and that together with 64 ( 04 ) it increases Acknowledgments 

with increasing magnetic field. 


IV. SUMMARY AND CONCLUSIONS 


In this paper we studied QCD FESR for the axial- 
vector current correlator in the presence of a magnetic 
field in the weak field limit eB < Sq, where Sq is the 
squared energy threshold for the onset of pQCD. Given 
that we work in the massless quark limit, the weak field 
limit can be understood as an expansion in the small pa¬ 
rameter eB/ sq. We have shown that the presence of the 
field modifies both the pQCD as well as the hadronic sec¬ 
tors of the FESR. The direct modification to the hadronic 
contribution for the equation that governs the behavior 
of So vanishes non-trivially. The only change in this sec¬ 
tor comes from the dependence of the pion decay con¬ 
stant on the magnetic field. In turn, this pion decay con¬ 
stant is proportional to the quark condensate through the 
GMOR relation. We have shown that in the limit where 
the magnetic field is small compared to the physical pion 
mass squared, the proportionality constant between 
and (qq) can be approximated to be the same as in vac¬ 
uum. Use of lattice QCD data on the quark condensate 
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Appendix A: Imaginary part of the hadronic 
contribution to the axial-vector current correlator 


The charged axial-vector current correlator in the pres¬ 
ence of a constant magnetic field, and in the limit q± —> 0 
is written as 




T) = -27E — 


<-i)‘ 


1=0 


s — (21 + l)eB 


(50) 
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The above expression can be split into two sums, one for 
the even and the other for the odd values of l, namely 


n ““M - - 2 4 E 

I l— 0 ,even 


1 )eB 


E- 


l—odd 


s-(2l + 1 )eB J ' 


(51) 


Therefore, by restricting the analysis to the region eB < 
so < 3 eB we can compute the discontinuity, or imaginary 
part of Eq. (15T1) , with the result 

Imn» AD (s) = /><5(s- eJ B), (55) 


where since s is strictly larger or equal to 0 , one has 


lim 7 -—-■„ — 

«->o (s — eB) + e 2 2 


-S(s-eB). 


(56) 


Pulling out a factor —1/AeB from both sums and adding 
and substracting the element with / = —1 we obtain from 
Eq. (El) 


n o AD (s) 


= Vi 


OO 1 

V_- 

' u _ 


AeB z —' // 


/eB -1 
4 


1 °° 1 

_ V_- 

eB ^ l' iL 


AeB ' /' _ s / e - H+1 s + eB 

l'= o L 4 . 


(52) 


where we defined V = 4 f° r the sum with even l and 
/' = hhl for the sum with odd l. The sums in Eq. m 
are divergent. In order to extract the finite piece we 
regularize them as 


n- D (s) 


2 fi lim 
1 


OO 

V_-_ 

' fj, _ s/eB-l y 


1 1 4eB (*' - 


V — 

^ (V 


l '=0 
1 


1 


4ei? ' d' _ s / e5+1 'je 
v =o \ L 4 ) 


s + eB 


= Vl lim 


e->i 1 AeB 


C(e, (s/eB - l)/4) 


AeB 


C(e, (s/eB + l)/4) - 


1 


s T eB 


(53) 


Finally, in the limit eB —> 0 the imaginary part of the 
correlator becomes 

Imllg AD (s) = />(*), (57) 

which coincides with the known value in the absence of 
a magnetic field. 


Appendix B: First order magnetic correction to the 
pQCD contribution to the axial-vector current 
correlator 


The contribution to the pQCD axial-vector current 
correlator of order (eB) is given by 


where 

n<i"V) 


n«'V) 


nju? 2 ) = + n^V) (58) 


= N c (q u B) 
= N c (q d B) 


d 4 k Tr [7r7 2 [7 ■ (k - q)\\h/i^u] 
(2 tt) 4 k 2 (k~q) 4 

d 4 k Tr [(# - ^ 7^7172 [7 • g]||7^] 
(27 t) 4 (k — q) 2 k 4 

(59) 


where £(a,z) is the Hurwitz zeta function. Expanding 
around e = 1 we find 


n o A °( s ) = Vi 


1 


1 


AeB 

1 


AeB 

1 


1 


e- 1 


e- 1 




-i> 


— ( s/eB — 1) 


-(s/eB + 1) 


s + eB 
2 


= -vi 

i 


4eB 


ip 


— (s/eB — 1) 


AeB 


tp 


-(s/eB + 1) 


eB 


(54) 


Note that the traces in Eqs. (USD are equal except for an 
overall sign, due to the ordering of elements inside the 
trace. Therefore, we evaluate one of the traces, 


Tr [7172 [7 -(k-q) || ]7m #7'«/] = k a (k - q) fj* 


= Tr 




= -ik a Tr [7 M 7a7^7p75 
X {& p [(fc — q) • it] 

- u p [(k - q) ■ b}}}, (60) 


where u p and b p are four-vectors describing the parti¬ 
cle’s rest frame and the direction of the magnetic field, 
respectively. In the rest frame, there are given by 


where ip(x) is the di-gamma function. We note that the 
divergent pieces cancel when e —> 1. Recall that ip(x) is 
singular for x = 0, —1, —2,.... In the region 0 < eB < sq, 
neither of the di-gamma functions in Eq. (1541) becomes 
singular. The first singularity for if>(—(s/eB — l)/4) hap¬ 
pens at s = eB and for ip(—(s/eB + 1)/4) at s = 3 B. 


Using that 


u = ( 1 , 0 ,0,0) 


b = (0,0,0,1). 

(61) 

Tr [ 7 ^ 7 ^ 75 ] = -Aie^vp, 

(62) 
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we can write Eq. as 

Tr [ 7172(7 ' (k - q) ||b/x#7*/] = 

-4fc% aj/p {b p [(k - q)-u\- u p [{k - q) ■ b}} . (63) 


Notice that when Eq. (l63l) is contracted with q fJ, q l/ the 
coefficient of the longitudinal structure vanishes. 
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